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Abstract 

We study the charged current Drell-Yan process and we evaluate the proton parton densi¬ 
ties uncertainties on the lepton transverse momentum distribution and their impact on the de¬ 
termination of the W boson mass. We consider the global PDF sets CTIO, MSTW2008CPdeut, 
NNPDF2.3, NNPDF3.0, MMHT2014, and apply the PDF4LHC recipe to combine the individual 
results, obtaining an uncertainty on mw that ranges between ±18 and ±24 MeV, depending 
on the final state, collider energy and kind. We discuss the dependence of the uncertainty 
on the acceptance cuts and the role of the individual parton densities in the final result. We 
remark that some PDF sets predict an uncertainty on m^ of ©(lOMeV); this encouraging 
result is spoiled, in the combined analysis of the different sets, by an important spread of the 
central values predicted by each group. 
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1 Introduction 


The very accurate measurement of the W boson mass performed at the Tevatron experiments 
CDF (m^ = 80.387 ± 0.019 GeV) and DO = 80.375 ± 0.023 GeV) [^, with a world 
average now equal to rriw = 80.385 ± 0.015 GeV [^, offers the possibility of a high-precision test 
of the gauge sector of the Standard Model (SM). There are prospects of a further reduction of 
the total experimental uncertainty, with a hnal error of 0(10) MeV for the combination of LHG, 
Tevatron and LEP results [^|^. 

The current best prediction in the SM is = 80.357 ± 0.009 ± 0.003 GeV and has been 
computed including the full 2-loop corrections [^, augmented by higher-order QGD corrections 
and by resumming reducible contributions. The uncertainty on this evaluation is mostly due 
to parametric uncertainties of the inputs of the calculation, the top mass value, the hadronic 
contribution to the running of the electromagnetic coupling, and to theoretical uncertainties. 

The simultaneous indirect determination of the top quark mass and the W mass, together with 
the direct determination of the Higgs boson mass, provides an important consistency check for the 
Standard Model [^; in turn the comparison of an accurate experimental Mw measurement with 
the predictions of different models might provide an indirect signal of physics beyond the SM. 

The W boson mass is extracted by means of a template £t technique applied to different 
observables of the charged-current (GG) Drell-Yan (DY) process, namely the lepton and neutrino 
transverse-momenta (p^, p^) and the lepton pair transverse mass rriT (see, for instance, 0|)' 
The differential distributions are computed with Monte Garlo simulation codes for different values 
of mw and are subsequently compared with the corresponding data: the value which maximizes 
the agreement is chosen as preferred value for mw- 

The present GDF and DO results are affected by a systematic error obtained as the combination 
of several elements, both of experimental and theoretical origin (see Tables IX and X in ref. 
for GDF and Table VI in ref. for DO). 

Among the experimental items, the most problematic ones are the determination of the lepton 
energy scale (0(7—17) MeV) and of the recoil scale and resolution (0(5—8) MeV), where the recoil 
ut is defined as the sum of the momenta of all the measured charged tracks, with the exception of 
the ones associated to the lepton(s). The mismeasurement of the recoil affects the determination 
of the W transverse momentum and the application of the selection cuts on this variable; in turn 
it affects the determination of the leptons transverse momenta, with an impact on the hnal mw 
value. The largest contribution of theoretical systematic error is due to the parametrization of 
the proton parton density functions (PDFs), which will be the main subject of the present paper, 
while another theoretical item present in these tables is the size of the missing QED effects not 
included in the available simulation tools, estimated to be of 0(4) MeV. 

The recoil modeling is an important element of the analysis of charged-current DY events, for 
it enters in the determination of the neutrino and of the charged lepton transverse momenta. The 
model is validated on neutral-current DY data of the lepton-pair transverse momentum distribu¬ 
tion, thanks to the fact that, in this latter case, the full information about the kinematics in the 
transverse plane can be reconstructed. The propagation of this calibration of the recoil model to 
the mw measurement is estimated to be of 0{3) MeV. The description of the recoil was optimized 
at the Tevatron experiments in the region of small lepton-pair transverse momenta, so that the 
DY events used for the mw determination are eventually selected imposing a cut ut < 15 GeV. 
In the following we will discuss the impact of this selection criterium, which will implemented in 
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our analysis as a cut on the transverse momentum of the lepton-neutrino pair. 

The theoretical contribution to the hnal systematic error enters in the analysis of the data via 
the templates. In fact the latter are computed with codes based on perturbative calculations, which 
are truncated at a hnite order in the expansion parameter. The proton PDFs, which are used to 
describe the partonic content of the proton, are affected by the error of the data from which they are 
extracted; this error propagates to the templates, inducing an additional systematic uncertainty 
of the htting tool. The interplay between the W and Z transverse momentum distributions and 
the proton PDFs has been discussed in ^ 12 , without a discussion on the consequence for the 
m^y determination. 

The aim of the present paper is to provide a quantitative assessment of the error induced by 
our imperfect knowledge of the proton PDFs, in the preparation of the templates used to £t niy,, 
in the case of the charged-lepton transverse momentum (p(p) distribution, as it is measured in the 
hadron collider processes pp —)■ lu + X. 

In 13 an analysis based on the lepton-pair transverse mass shows that PDF uncertainties do 


not challenge a measurement of at 10 MeV accuracy. As discussed above, the measurement 
of the lepton pt has different, and to a certain extent complementary, systematic uncertainties, 
compared to the lepton-pair transverse mass. On the other hand, getting theoretical predictions 
for the pip distribution can be quite challenging due to the high sensitivity of this exclusive ob¬ 
servable to the details of the description of the radiation emitted. In this respect, also the PDFs 
parametrization has a direct impact on the shape of the distribution and, in turn, on rriyy. 


1.1 The lepton transverse momentum distribution and m 


w 


The sensitivity of the lepton transverse momentum distribution to the precise value of the W 
boson mass is due to its Jacobian peak which has its maximum for ~ m^/2. A change of rriw 
in the simulation codes by 2,10, 20 MeV, with respect to a fixed reference value, yields a distortion 
of the distribution in the per mill range, as illustrated in Figure Any effect (perturbative QCD 
corrections, PDF uncertainties, etc.) that induces a change of the shape of this distribution of 
similar size represents a source of systematic theoretical uncertainty on the determination; in 
particular a measurement at the 10 MeV level requires a control of the shape of the templates at 
the 1 per mill level or better. 

There are two mechanisms that yield a distribution of the lepton transverse momentum in the 
DY processes: the decay of the gauge boson and its recoil against QCD (and in smaller amount 
QED) radiation. The initial state radiation collinear divergences make any hxed-order prediction 
for this quantity unreliable, because of the important contributions in the region of small gauge 
boson transverse momenta, which have to be resummed to all orders. It is thus necessary to 
use a code that implements the resummation to all orders of multiple gluon emissions, either 
analytically or in a numerical approach via Parton Shower (PS), to obtain a physically sensible 
prediction. In this study we use the POWHEG Monte Carlo event generator 21 , matched with 


the PYTHIA 22 QCD PS. The accuracy of this code on the inclusive DY cross section is NLO- 


QCD, while from the point of view of the enhancement due to the logarithms of the transverse 
momentum of the lepton pair, the lepton transverse momentum distribution has LL accuracy. We 
consistently choose NLO-QCD PDF distributions. The POWHEG event generator is currently used 
by the ATLAS and CMS (see e.g. respectively refs. and [^) collaborations to study the DY 
processes and provides a good description of the data. We thus consider it as a valid starting point 
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Figure 1: Ratio of lepton transverse momentum distributions which have been generated with 
different W boson masses. 


to study the propagation of the PDF uncertainty in the preparation of the templates eventually 
used to £t the data. 

A more accurate description of the charged-lepton transverse momentum distribution can be 
obtained by the use of codes that include higher-order QCD corrections, like ResBos 26,27 or 
DYRes 25 , or EW and QED multiple photon effects, like POWHEG 28-^. As it is well known, 
hnal state QED radiation plays a crucial role in the precise determination of MW 31 . However, 
since in this study our main focus is on the assessment of the impact on the m^y determination of 
the PDF uncertainty, we choose a fast code that yields a basic realistic description of the shape 
of this distribution. The difference with the predictions that one could obtain adopting one of 
the other above listed codes belongs to the class of mixed PDF x higher-order effects and it can 
be estimated as a perturbative correction to the results of the present study; the use of a code 
that includes hnal state QED effects modihes the basic shape of the templates, with a shift of the 
central rriw value of 0(200) MeV (^; our ansatz, in the absence of an explicit check, is that this 
modification of the shape yields also a rescaling of all the PDF uncertainties at most of 0(10%) 
of the results discussed in this paper; this estimate follows from the comparison between the QED 
shift and the size of the region of the distribution sensitive to an m^y variation, which is of at 
least of one, but more realistically of a few W decay widths. We would thus obtain a change in our 
results, for a given PDF set, by 0(2 MeV), which is comparable to the statistical accuracy that 
we can claim in the template £t and to the error, of experimental origin, on the PDF uncertainty 


itself 33 


The plan of the paper is the following: in Section we recall the dehnition of some basic 
theoretical tools which will be used in the study; in Section we present our numerical results, 
discussing the PDF uncertainty on the lepton transverse-momentum distribution and on the myy 
determination; we also consider the dependence of the niw PDF uncertainty on the acceptance 
cuts and comment on possible future developments; in Section we draw our conclusions. 
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2 Theoretical tools 


In this section, we briefly ontline the strategy adopted to estimate the PDF nncertainty in the 
determination of mw at hadron colliders: we refer the interested reader to 13 for more details. 


2.1 Template fit 

In this paper we discnss the nncertainty on rriyy of different PDF sets. In order to make a 
qnantitative evalnation, we follow some basic steps: 


1. we generate the lepton transverse momentnm O with different PDF replicas, keeping the W 
mass fixed at a given valne and we treat each distribntion as a set of pseudodata] 


2. we compnte the templates, i.e. the distribntions that are nsed to fit the psendodata, with 
one specific choice for the PDF set (NNPDF2.3, replica 0) and we let rriyy assnme all the 
valnes in the interval [80.312,80.470] in steps of 2 MeV; 


3. we compare a given psendodata distribntion, obtained with a given PDF replica labelled by 
i, with all the templates labelled by j] in each comparison we compnte an indicator 



1 


iV, 


bins 


-^bins 


k=l 


template 


^t,aata \ 

k J 


^^i^datay _|_ (^^Mempiatey _ 2Cov{data, template) 


( 1 ) 


where Ok and ak are respectively the valne of the distribntion and its associated error in the 
bin k and Cov is the covariance between the two distribntions; 

4. the template j that yields the minimnm valne of Xij is the one that best describes the 
psendodata and its associated m^j valne is thus the preferred value associated to the replica 
i] the difference Amw,i = — m^o, is the shift induced by the PDF replica i chosen for 

that set of psendodata; in other words it is the difference between the results that we would 
obtain when fitting the real data if we prepared the templates with the replica i instead of 
the replica 0 of NNPDF2.3. 


2.2 PDF uncertainties 

The proton PDF sets considered in this study are MSTW2008CPdeut 


14 , CTIO 15 , NNPDF2.3 16 


NNPDF3.0 17 and MMHT2014 18 . They are called global sets because they include all the available 


relevant hard scattering data. Each collaboration provides a prescription to estimate the PDF 
uncertainties!^ in particular we recall the formula for the symmetric error in the Hessian approach 
(CTIO, MSTW2008CPdeut, MMHT2014) for a generic observable X 


AX = 


N, 


eigenvectors 


\ h 


E Pt-AT]' 


( 2 ) 


'^We refer to the original publications for more details. 
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where the sum runs over the Neigenvectors eigenvectors in parameter space, with the associated 
pairs of replicas (+ and -). Instead with NNPDF the average and the standard deviation over the 
ensemble {q} of Nj-ep PDF replicas provide the estimate of the best value and of the error on the 
observable J^: 

Nrep 

= • (3) 

1\ r - - 


k=i 


(yj = 


iv, 


vep 


Nrep 


k=l 


1/2 


(n{9'‘>}] - (n{9}]))' 


(4) 


The results obtained with these three PDF sets can be combined according to the current PDF4LHC 


recommendation 20 , to hnd a conservative estimate of the PDF uncertainty. 


In this paper we apply this procedure to two observables, namely the lepton transverse mo¬ 
mentum distribution and the W mass determined with the template £t procedure. 

2.3 Correlation functions 

A useful quantity to evaluate the role of the different parton densities in the hadronic cross 
section is the correlation function p between the parton-parton luminosities and the charged- 
lepton distribution at a given value of the transverse momentum. The parton-parton luminosity 
is dehned as Vij{x,T) = fi{x, pp) where fi{x,p?p) is the density describing a parton i at 

a scale pp and t = ^ with M the final state invariant mass and S the hadronic Mandelstam 
invariant. The correlation p is dehned as 


p{x,t) = 




(J- 


PDF^PDF 


V 


^3 


(P \ 


(5) 


da / 

where the angle brackets indicate average with respect to the different PDF replicas. 
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3 Numerical results 


3.1 Input parameters and setup 


We simulate the processes pp —)■ W~^ —)■ + X and pp —)■ W~ —j- + X in proton- 

antiproton collisions with -s/S = 1.96 TeV and in proton-proton collisions with \/~S = 8,13, 33,100 
TeV energies. In the absence of QED effects, not considered here, our results will be identical to 
those obtained with electrons instead of muons. We consider the PDF sets MSTW2008CPdeut 14 , 
CTIO [^, NNPDF2.3 [^, NNPDF3.0 17, MMHT2014 18, and use the corresponding values of 

as{m, 


(-) 


We use the following values for the input parameters in the Monte Carlo codes: 


= 1.16637 10-5 GeV-2 

mw = 

-- 80.398 GeV 

rriz 

= 91.1876 GeV 

sin^ 9w = 1 — m?^lm\ 

Tw = 

2.141 GeV 

Tz 

= 2.4952 GeV 

Vai = 0.222 

Vcs = 

0.975 

Vcb 

= 0 

Kd = 0.975 

Vus = 

0.222 

Vub 

= 0 

Vtd = ^ 

Vts = 

0 

Vtb 

= 1 


The charm quark in the partonic cross section is treated as a massless particle, while the bottom 
quark does not contribute because of the vanishing top density in the proton. As for the kinematic 
cuts, we used those summarized in Table [T| similar to those used in the corresponding experimental 
analysis: the main difference between the Tevatron and LHC is the wider acceptance for the 
rapidity of the leptons in the latter case. The p^j, distribution has been studied in the interval 29 
GeV< p\ < 49 GeV, with a bin size of 0.5 GeV. All the following analyses are performed with 
bare leptons both in the pseudodata and in the templates. 


Tevatron 

LHG 

> 25 GeV 
$T > 25 GeV 
\9^ < 1-0 
pf < 15 GeV 

p^ > 25 GeV 
Ipj' > 25 GeV 
\r|^ < 2.5 
p^ < 15 GeV 


Table 1: Selection criteria for DY W ^ lu events for the Tevatron and the LHC. 


The Monte Garlo simulation requires a specihc, technical comment. The effects under study 
are deformations of the shape of the lepton transverse momentum distribution at the per mill level, 
either due to a variation of the mw value or to a different PDF replica choice. This distribution 
receives contributions from a large fraction of the available hnal-state phase space, making very 
difficult an accurate dedicated sampling. As a consequence, Monte Garlo statistical fluctuations 
at the per mill level are present also with hundreds of milions of simulated unweighted events. 
The solution to this problem is found using a reweighting technique, based on the remark that 
both the dependence on the PDFs and the dependence on mw factorize from the rest of the hxed- 
order partonic cross section. Only one simulation, one sequence of events is used to generate all 
the templates and all the pseudodata: the weight wq associated to each event is corrected by an 


7 








appropriate reweighting factor to account for different replica or, separately, mw value choices. 


Wq —)■ Wj = Wq- 


Wq -)■ Wi = Wq 


(s — m 


wo) 




fi{xi)gi{x2) 


2 

W,j 


Jo 


ixi)9b 


ix2) 


template j 
replica i 


( 6 ) 


where /, g are two generic parton densities. In the POWHEG formulation this rescaling spoils the 
exact NLO accuracy of the hnal result, by terms generated by the POWHEG Sudakov form factor. 
The size of the latter could not be distinguished from Monte Carlo statistical fluctuations, when 
we compared two distributions, one obtained with an exact simulation and the other introducing 
the new PDF replica via the above rescaling. The main results of this study should thus not be 


affected. An update of the POWHEG generator is in progress 19 to restore the NLO accuracy after 
reweighting. 

Since the events used are exactly the same, the statistical fluctuations of the different dis¬ 
tributions (templates and pseudodata) are highly correlated (correlation is about 0.987, almost 
constant over the bins, with SOM of events) and cancel to a large extent when we compute the 
difference of two cross sections in a bin. A statistical error of ±2 MeV can be obtained with the 
simulation of 1620 millions of events. We checked that, with increasing statistics, the result of 
each individual ht is stable, because the uniform reduction of the statistical error in the different 
p(j_ bins. We stress that, rather than an absolute prediction of the mw value, we are interested in 
the quantitative assessment: i) of the relative difference between various PDF sets and ii) of the 
PDF uncertainty within one set, which is dehned as the spread with respect to a central value. In 
both cases we have to provide a solid estimate of differences bin by bin and the reweighting pro¬ 
cedure allows to efficiently remove the Monte Carlo fluctuation effects, leaving only the physically 
relevant shifts. 


3.2 PDF uncertainty of the distribution 

We study the percentage PDF uncertainty on the lepton transverse momentum distribution and 
also on the associated normalized distribution dehned as 


da 1 da 



(7) 


As it is well known from the Tevatron experiments 0i. the uncertainty is in general reduced 
in the normalized observable. In fact, each PDF replica in a set contributes in a different way 
to the shape and to the overall normalization of the physical distributions; by considering the 
normalized distributions of Equation we are sensitive mostly to the shape change; the latter 
is the most relevant item in the determination of mw, because we associate the precise position 
and shape of the Jacobian peak to the value of the gauge boson mass. In Figure we show the 
PDF uncertainty of the lepton transverse momentum distribution and of the associated normalized 
distribution, computed at the Tevatron, at the LHC 8 and 13 TeV with different PDF sets, both in 
the cases of W~^ and W~ production, in presence of the additional cut < 15 GeV with respect 
to the choices indicated in Table The percentage uncertainty of the normalized distributions is 
at the few per mill level at the Jacobian peak, and could mimic the effect of a mw shift by (P(IO) 
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Figure 2: Percentage size of the PDF uncertainty on the lepton transverse momentum distribu¬ 
tion, computed with different PDF sets. In addition to the basic acceptance criteria of Table [1} 
a cut pY < 15 GeV on the lepton pair has been applied. The lower lines refer to the normalized 
distributions of Equation the upper lines to the standard ones. Results for W~^ (left) and W~ 
(right) production at the LHC 8 TeV (middle plots) and 13 TeV (lower plots); results for the 
Tevatron in the upper plot. 


MeV. In Figure we use the PDF set NNPDF3.0 and study the change of the PDF uncertainty 
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Figure 3: Percentage size of the PDF uncertainty on the lepton transverse momentum dis¬ 
tribution, computed with the NNPDF3.0 set at the LHC at different energies. Results for W~^ 
(left) and W~ (right) production, in the case of absolute (upper lines) or normalized (lower lines) 
distributions. 

with the collider energy, in presence of the additional cut pY < 15 GeV. The uncertainty of the 
distribution increases, as function of the collider energy, from 1.5 to 2.5%, while in the normalized 
case the uncertainty is almost independent of the energy. 


3.3 Impact of the PDF uncertainty on the rriw determination 


The template £t procedure, described in Section |2.1[ has been applied to the distributions com¬ 
puted with all the replicas of the different PDF sets under study; the corresponding preferred rriw 


values have been combined, according to the rules described in Section 2.2, to derive the uncer¬ 


tainty on the niw extraction due to the PDFs. The £t interval has been chosen to be pj_ G [29,49] 
GeV, in order to minimize the contribution to the PDF uncertainty from the tails of the distri¬ 
bution above and below the Jacobian peak. The template fit has been applied to our pseudodata 
generated with a fixed value of the W boson decay width F^. We have checked that our results 
are weakly dependent on the choice of this parameter: we repeated the £t using for F^ a value 
modihed by iur, where ur = 0.042 GeV is the current experimental error, and we found that the 
prediction for the PDF uncertainty on rriw gets modihed by 1-2 MeV, depending on the selection 
cuts. The results for the Tevatron and for the LHG 8 and 13 TeV are presented in Table and 
are also summarized in Figures In the upper half of Table (and in Figure left plot) no 
additional cut on has been imposed on the lepton pair, whereas in the lower half of the same 
Table (and in Figure [fright plot) a cut p^ < 15 GeV has been applied. 

The PDF uncertainty rehects the experimental error of the data from which the parton densities 
are extracted and also the different methodologies used in their determination. As it can be 
observed from Figure the estimate of the PDF uncertainty on predicted by the different 
PDF collaborations differs by a factor up to 3 between the different groups. The uncertainty on 
extracted from normalized distributions, with the basic selection criteria of Table ranges 
from 12 to 23 MeV at the Tevatron, from 12 to 29 MeV at the LHG 8 TeV and from 11 to 34 
MeV at the LHG 13 TeV. Imposing on the lepton pair a cut p^ < 15 GeV modifies these results; 
the ranges of the PDF uncertainties become: from 11 to 17 MeV at the Tevatron, from 7 to 17 
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Figure 4: Summary of the PDF uncertainty on ruw computed with different PDF sets, colliders 
and final states. The basic acceptance criteria have been used in the left plot, while in the right 
plot an additional cut <15 GeV has been applied. 
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MeV at the LHC 8 TeV and from 6 to 18 MeV at the LHC 13 TeV. 

In addition, the PDF sets nnder stndy differ in the parametrization that they adopt to describe 
the proton strnctnre; the latter affects the best description of the PDFs and in tnrn the best 
prediction of the central rriw valne. The spread Agets of the central valnes, dehned as the difference 
between the largest and the smallest central valnes, is a second component of the hnal PDF 
nncertainty on 

A conservative estimate of the nncertainty on m^r, that combines the two elements of nncer¬ 
tainty described above, can be obtained by compnting the envelope of the predictions nnder stndy, 
according to the PDF4LHC recipe and by measnring the half-width oi the resnlting band. 

We inclnde, in the evalnation of the envelope, the resnlts of the sets CTIO, MSTW2008CPdeut and 
NNPDF2.3, becanse they are based on the same sets of data, making their comparison homoge- 
neons. These resnlts are presented in Table We observe that the spread Ag^ts represents a large 
contribntion, np to 35% of the overall nncertainty . In Table we compnte the envelope of the 
resnlts obtained with two more modern PDF sets, namely NNPDF3.0 and MMHT2014, which inclnde 
recent pnblic data from the LHC. We observe that the width of the envelope ranges between 16 
and 32 MeV, depending on the collider energy and kind and on the hnal state; more interesting, 
the spread of the two central valnes is below 5 MeV in the W~ case at the LHC, while it is above 
15 MeV in the case and at the Tevatron. 

From Table we can appreciate the impact of the inclnsion of the new LHC data, which have 
been nsed in the determination of the NNPDF3.0 set. Beside a few MeV offset for the central 
valnes, it is possible to observe a small (few MeV) redaction of the PDF nncertainty, which is 
ronghly 20% smaller than the one compnted with NNPDF2.3. For MMHT2014 the nncertainties are 
similar or slightly larger than the ones obtained with MSTW2008CPdeut. 

The dependence of the PDF nncertainty with the collider energy is illnstrated in Table nsing 
the NNPDF3.0 PDF set. 
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no pY cut 

p'f < 15 GeV 


SpDF (MeV) 

(MeV) 

6pDF (MeV) 

(MeV) 

Tevatron 1.96 TeV 

27 

16 

21 

15 

LHC 8 TeV W+ 

33 

26 

24 

18 

w- 

29 

16 

18 

8 

LHC 13 TeV W+ 

34 

22 

20 

14 

w- 

34 

24 

18 
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Table 2: Half-width 6pdf of the envelope of the PDF uncertainty intervals by CTIO, 

MSTW2008CPdeut and NNPDF2.3. Corresponding spread ^sets of the central predictions. 



no cut 

< 15 GeV 


6pdf (MeV) 

A,et. (MeV) 

SpDF (MeV) 

A,ei. (MeV) 

Tevatron 1.96 TeV 

16 

4 

13 

9 

LHC 8 TeV W+ 

32 

33 

21 

21 

w- 

22 

6 

12 

0 

LHC 13 TeV W+ 

30 

24 

18 

16 

w- 

23 

16 

11 

5 


Table 3: Same as in Table now considering only the two recent PDF sets NNPDF3 .0 and 
MMHT2014. 


normalized distribution, additional cut pY < 15 GeV 


8 TeV 

13 TeV 

33 TeV 

100 TeV 

W+ 

80.395 ± 0.009 

80.400 ±0.010 

80.402 ±0.010 

80.404 ± 0.013 

W- 

80.398 ± 0.007 

80.391 ±0.006 

80.385 ±0.007 

80.398 ±0.011 


Table 4: Estimate of the central values and of the PDF uncertainty on rriyy, extracted from the 
lepton transverse momentum distributions simulated with the NNPDF3.0 set at different proton- 
proton collider energies. 


12 
































absolute distributions 

collider / channel 

CTIO 

MSTW2008CPdeut 

NNPDF2.3 

NNPDF3.0 

MMHT2014 

Tevatron, 

80.406 + 0.043 - 0.046 

80.428 + 0.025 - 0.017 

80.400 + 0.030 

80.427 + 0.018 

80.430 + 0.022 - 0.022 

LHC 8 TeV, W+ 

80.394 + 0.040 - 0.029 

80.422 + 0.025 - 0.016 

80.398 + 0.020 

80.406 + 0.019 

80.428 + 0.027- 0.022 

w- 

80.444 + 0.055 - 0.062 

80.390 + 0.038 - 0.036 

80.398 + 0.030 

80.441 + 0.027 

80.404 + 0.041 - 0.048 

LHC 13 TeV,H"+ 

80.396 + 0.045 - 0.034 

80.416 + 0.020 - 0.020 

80.398 + 0.022 

80.414 + 0.022 

80.422 + 0.030 - 0.024 

w- 

80.416 + 0.088 - 0.065 

80.374 + 0.044 - 0.033 

80.398 + 0.031 

80.426 + 0.037 

80.384 + 0.037- 0.049 

normalized distributions 

collider / channel 

CTIO 

MSTW2008CPdeut 

NNPDF2.3 

NNPDF3.0 

MMHT2014 

Tevatron, 1T+ 

80.400 + 0.022 - 0.025 

80.414 + 0.016-0.016 

80.398 + 0.012 

80.408 + 0.013 

80.412 + 0.014-0.010 

LHC 8 TeV, W+ 

80.398 + 0.032 - 0.026 

80.424 + 0.014- 0.019 

80.398 + 0.016 

80.395 + 0.014 

80.428 + 0.016 - 0.024 

w- 

80.416 + 0.026 - 0.025 

80.398 + 0.011 - 0.014 

80.398 + 0.014 

80.396 + 0.012 

80.402 + 0.019 - 0.024 

LHC 13 TeV,H"+ 

80.406 + 0.039 - 0.029 

80.420 + 0.017- 0.014 

80.398 + 0.018 

80.404 + 0.016 

80.428 + 0.020 - 0.026 

w- 

80.422 + 0.030 - 0.023 

80.398 + 0.008 - 0.015 

80.398 + 0.015 

80.386 + 0.011 

80.402 + 0.019 - 0.024 

absolute distributions, additional cut < 15 GeV 

collider / channel 

CTIO 

MSTW2008CPdeut 

NNPDF2.3 

NNPDF3.0 

MMHT2014 

Tevatron, W~^ 

80.412 + 0.024 - 0.024 

80.424 + 0.018 - 0.017 

80.398 + 0.013 

80.420 + 0.014 

80.426 + 0.009 - 0.021 

LHC 8 TeV, W+ 

80.392 + 0.026 - 0.021 

80.414 + 0.020 - 0.011 

80.398 + 0.014 

80.403 + 0.014 

80.418 + 0.019-0.017 

w- 

80.422 + 0.039 - 0.034 

80.394 + 0.019 - 0.023 

80.398 + 0.017 

80.423 + 0.017 

80.400 + 0.023 - 0.028 

LHC 13 TeV,H"+ 

80.392 + 0.028 - 0.022 

80.410 + 0.012-0.016 

80.398 + 0.014 

80.408 + 0.014 

80.414 + 0.016-0.019 

w- 

80.408 + 0.042 - 0.037 

80.386 + 0.019 - 0.021 

80.398 + 0.016 

80.410 + 0.018 

80.388 + 0.021 - 0.025 

normalized distributions, additional cut < 15 GeV 

collider / channel 

CTIO 

MSTW2008CPdeut 

NNPDF2.3 

NNPDF3.0 

MMHT2014 

Tevatron, W~^ 

80.400 + 0.018 - 0.016 

80.414 + 0.013-0.015 

80.398 + 0.009 

80.403 + 0.011 

80.412 + 0.006 - 0.012 

LHC 8 TeV, W+ 

80.396 + 0.017- 0.018 

80.414 + 0.012-0.011 

80.398 + 0.010 

80.395 + 0.009 

80.416 + 0.011-0.014 

w- 

80.406 + 0.016 - 0.011 

80.398 + 0.005 - 0.012 

80.398 + 0.010 

80.398 + 0.007 

80.398 + 0.008 - 0.016 

LHC 13 TeV,H"+ 

80.400 + 0.020 - 0.017 

80.412 + 0.010-0.011 

80.398 + 0.011 

80.400 + 0.010 

80.416 + 0.010-0.015 

w- 

80.408 + 0.017- 0.009 

80.396 + 0.010 - 0.006 

80.398 + 0.009 

80.391 + 0.006 

80.396 + 0.009 - 0.013 


Table 5: Estimate of the central values and of the PDF uncertainty on mw, extracted from the lepton transverse momentum 
distributions simulated with different PDF sets and acceptance cuts. The templates have been generated with NNPDF2.3 replica 
0. The pseudodata for the different PDF sets have been simulated by setting ttiw = 80.398 GeV. 




3.4 PDF uncertainty dependence on the acceptance cuts 

The results presented in Section |3.3| have been obtained imposing on the leptons the basic cuts 
of Table The dependence of the PDF uncertainty on additional cuts on the lepton-pair 
transverse momentum or on the charged-lepton pseudorapidity acceptance interval is presented 
in Table This study suggests possible optimizations of the event selection, to minimize the PDF 
uncertainty impact. It also offers a link to the dependence of the PDF uncertainty on the different 
flavors in the proton and on the most problematic range in partonic-x. 


normalized distributions 

cut on p^ 

cut on \pi\ 

GTIO 

NNPDF3.0 

inclusive 

h\ < 2.5 

80.400 + 0.032 - 0.027 

80.398 + 0.014 

p^ < 20 GeV 

h\ < 2.5 

80.396 + 0.027 - 0.020 

80.394 + 0.012 

p^ < 15 GeV 

h\ < 2.5 

80.396 + 0.017- 0.018 

80.395 + 0.009 

p^ < 10 GeV 

h\ < 2.5 

80.392 + 0.015 - 0.012 

80.394 + 0.007 

p 7 < 15 GeV 

\vi\ < 1.0 

80.400 + 0.032 - 0.021 

80.406 + 0.017 

p 7 < 15 GeV 

h\ < 2.5 

80.396 + 0.017- 0.018 

80.395 + 0.009 

p 7 < 15 GeV 

h\ < 4.9 

80.400 + 0.009 - 0.004 

80.401 + 0.003 

p 7 < 15 GeV 

1.0 < \pi\ < 2.5 

80.392 + 0.025 - 0.018 

80.388 + 0.012 


Table 6: LHC 8 TeV, production. Impact of different acceptance cuts. The two cuts > 25 
GeV and.^r > 25 GeV are always applied. In the hrst four rows we vary the cut on for hxed 
\rji\ interval. In the second four rows we vary the pseudorapidity acceptance, with < 15 GeV. 



Figure 5: Shape of the normalized differential distribution da/dx for different p^ cuts (left plot). 
Ratio of the previous shapes with different p^ cuts with respect to the inclusive (no p^ cut) 
distribution (right plot). 


We observe that the region at large p^f yields an important contribution to the PDF uncer¬ 
tainty, which can be reduced by a suitable cut on this variable. A tight cut like p^ < 10 GeV 
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Figure 6: Shape of the lepton transverse-momentum (left panel) and of the lepton pseudorapidity 
(right panel) distributions, in presence of different additional cuts on the lepton-pair transverse 
momentum ■ 



Figure 7: Percentage uncertainty of the individual parton densities f{x,m‘^) of NNPDF3.0 (left 
plot). Correlation of different parton-parton luminosities with the charged-lepton distribution 
at p(i_ = 40.5 GeV, computed with different acceptance cuts on \rii\ and with < 15 GeV. 


could bring the uncertainty below the 10 MeV level. The experimental problem to accurately 
select the events that pass the cut can be a limiting factor for the improvement in this direction. 

The impact of the cut on the lepton-pair transverse momentum can be explained by studying 
the change of the relative contribution of the medium- vs the large-x PDF region, where x is the 
fraction of momentum of the parent hadron carried by the incoming parton. In Figure]^ (left plot) 
we show the normalized da/dx distributions, where x is the fraction of longitudinal momentum 
carried by the partons of one given hadron in the scattering they are computed with different p'^ 
cuts, and express the relative contribution of a given partonic x to the cross section. In Figure 

^The choice of the hadron is not relevant, because the contribution of the partonic subprocesses is symmetric 
for exchange of hadrons 1 and 2 
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Figure 8: Percentage PDF uncertainty of the charged-lepton distribution (left plot) and shape 
of the differential distribution da/dx (right plot), computed with different acceptance cuts on |n;| 
and with pY < 15 GeV. 


(right plot) we show the ratio of the previous distributions, computed with different p^ cuts, with 
respect to the inclusive (no p^ cut) normalized distribution. These ratios express the relative 
change of the weight of the various x intervals, in presence of a cut. We thus recognize that the 
pY < 15 GeV cut enhances the x < 0.004 region and suppresses the contribution at a; > 0.004. 
Since the PDF uncertainty of all the densities rapidly increases for x > 0.1 (cfr. Figure]^ left 
plot), the effect of the pY cut is a reduction of the global PDF uncertainty affecting the rriw 
determination. A second effect of the cut is a change of the basic shape of the distribution, which 
becomes steeper and closer the LO one, above the Jacobian peak, as it is shown in Figure this 
modihcation increases the sensitivity of the htting procedure, which becomes more stable, because 
large shifts are more penalized with respect to the case of a broader distribution. In right panel of 
Figure]^ we show the normalized lepton pseudorapidity distribution, computed for different values 
of the pY cut. We observe that with tighter cuts the distribution develops two peaks at forward 
and backward rapidities. These regions are dominated by the contribution of at least one valence 
quark, whose PDF uncertainty is smaller than the one of the corresponding sea component. 

We observe that, for hxed cut on pY■, the PDF uncertainty decreases from 17 (26) to 3 (6) 
MeV with NNPDF3.0 (CTIO), as one enlarges the charged-lepton rapidity cut, from 1.0 to 4.9. This 
reduction is consistent with the smaller PDF uncertainty of the lepton transverse momentum 
distribution with the cut \rii\ < 4.9 shown in Figure]^ (left plot). In this case the problematic 
point is the possibility of an accurate measurement of the lepton properties in the large rapidity 
regions of the detector. 

The impact of the cut on the charged-lepton pseudorapidity can be explained hrst of all 
by recalling that a lepton transverse momentum distribution fully integrated over the lepton- 
pair rapidity (without acceptance cuts) would depend on the PDFs only via a single numerical 
factor, which drops out when we study the normalized distributions. This ideal limit can be 
reached, in a realistic setup, by enlarging the charged-lepton pseudorapidity acceptance. More in 
detail, with different maximal values of rji, we observe a corresponding change of the shape of the 
da/dx distribution, shown in Figure]^ (right plot): the bulk of the distribution is peaked around 


16 


























5 ■ 10“^, 1 ■ 10“^, 5 ■ 10“^ respectively for \rii\ < 1, 2.5, 4.9. 

First of all we observe in Table that the two PDF uncertainties on niw extracted imposing 
the cuts \r]i\ < 1.0 or 1.0 < \r]i\ < 2.5 are separately larger than the one obtained with \r]i\ < 2.5. 
Indeed, the sensitivity to partonic x obtained by varying the cut on rji makes evident the presence 
of the momentum sum rules, which have to be fulhlled by all the replicas: the more inclusive setup 
is thus more stable with respect to a PDF replica variation than the more exclusive cases. This 
uncertainty reduction is even more pronounced with \rii\ < 4.9. 

Second, in the region 5 ■ 10“^ ^ x < 1 ■ 10“^ the strange density has its maximal uncertainty, 
which is more than 3 times larger than the one of all the other parton densities, as shown in 
(left plot); in this same region the parton-parton luminosity cs has a weak positive correlation 
with respect to the PDFs, dehned in Equation with the charged-lepton transverse momentum 
distribution at = 40.5 GeV, so that its contribution to the cross section and to the total PDF 
uncertainty sums together with the ones of the other channels. In the interval 3T0“^ ^ x < 7T0“^ 
the strange density still has a PDF uncertainty 2.5 times larger than the others, but in this region 
the cs luminosity has a negative correlation with the distribution. In this case there are non trivial 
compensations between the contributions of the various partonic subprocesses, yielding a more 
stable result with respect to the PDF replica choice. This behavior of the parton densities, in the 
case of this quite inclusive observable, with the acceptance cut \r]i\ < 4.9, reflects the enforcement 
in the global £t of the sum rules that have to be satished by the PDFs. 


3.5 Comparison with previous studies 


The PDF uncertainty affecting the mw determination from the study of the lepton transverse 
momentum distribution has been estimated in and in to be respectively 12 and 11 MeV. This 
evaluation was based on the simulation code ResBos and on the use of the PDF sets cteqG . 6 
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and MSTW2008 We repeated the estimate of the uncertainty in the Tevatron setup also with 
CteqG. 6 , using POWHEG+PYTHIA and templates computed with NNPDF2.3 replica 0; we obtain 
mw = 80.396 -P 0.015 — 0.016, i.e. a slightly larger PDF uncertainty compared to the previous 
estimate. 

At variance with the transverse mass case 


13 , the study presented in this note, with events 


treated at generator level, is moderately sensitive to detector effects and should thus represent a 
realistic estimate of the overall size of the PDF uncertainty and of the relative behavior of the 
different PDF sets. 


3.6 The role of the lepton-pair transverse momentum distribution in 
NC-DY 

The description of the lepton transverse momentum distribution depends on the treatment of 
initial-state QCD radiation, to obtain a correct lepton-pair transverse-momentum distribution 
and in turn the correct contribution to the lepton transverse momentum. At low lepton-pair 
transverse momenta there are non-vanishing non-perturbative effects, which can be accounted 
for by means of ad hoc models, upon which the hnal result of mw depends. The uncertainties 
of the PDFs and of the modeling of an intrinsic component k± of the transverse momentum of 
the partons inside the proton are entangled in the lepton-pair transverse momentum distribution, 
because of the different contribution of the various flavors to the transverse momentum spectrum; 
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in other words, it is not possible, in principle, to derive a universal, flavor independent, model 
of the intrinsic k^. This statement has been investigated in the past (see e.g. [^) where the 
Tevatron data were described by a universal flavor-independent non-perturbative functions. 

A reduction of this dependence can be obtained by considering new observables, defined as 
ratios of the CC-DY observables with respect to their analogous ones in the case of NC-DY 36 


The similarities in the initial-state QCD radiation patterns determine a correlation between the 
CC-DY and the NC-DY quantities, which in turn yields a reduction of the error that affects 
the ratio. One should however keep in mind that it is not possible to expect a full correlation 
between the CC-DY and the NC-DY observables, because of the different flavor structure of the 
subprocesses in the two cases and because of the different phase-spaces available. 

Given the entanglement between PDF and intrinsic k±_ uncertainties, the estimate of the PDF 
uncertainty alone presented in this paper, for the CC-DY case, could be a slight overestimate of 
its contribution to the total non-perturbative uncertainty. On the other hand, the estimate of the 
PDF uncertainty for the ratio of CC-DY with respect to NC-DY observables should offer a more 
reliable result, thanks to the weaker model dependence. A detailed study of these ratios and of 
the associated theoretical uncertainties will be presented elsewhere. 


4 Conclusions 

We presented a quantitative assessment of the PDF uncertainty affecting the extraction of the 
W boson mass from the study of the charged-lepton transverse-momentum distribution in the 
charged-current Drell-Yan process, at different hadron colliders, for different collider energies. 
The study, conducted at generator level, is based on the Monte Carlo code POWHEG interfaced 
with the PYTHIA QCD parton shower and uses the NNPDF2.3 PDF set (replica 0) to prepare the 
templates used in the fitting procedure. The results are summarized in Figure]^ and in Table 
The study provides information about the relative distance between the NNPDF2.3 and the other 
sets considered (CTIO, MSTW2008CPdeut); this distance is expressed by the difference between the 
best predictions of the various sets and ranges between 8 and 15 MeV, depending on the collider, 
on the energy and on the final state considered; these results rely on the application of a cut 
on the lepton-pair transverse momentum, < 15 GeV. The study provides an estimate of the 
PDF uncertainty according to the prescriptions of each PDF group: the individual values range 
between 6 and 18 MeV, again depending on the considered setup and always in presence of the 
cut on the lepton-pair transverse momentum. The combination of the two previous uncertainties, 
according to the PDF4LHC recipe, leads to a global PDF uncertainty that ranges between 18 and 
24 MeV. The analysis of more modern sets, like NNPDF3.0 and MMHT2014, does not change this 
overall picture, but makes evident some differences in the description of with respect to W~ 
production. 

We remark that the differences between the PDF sets considered here are large compared to an 
accuracy goal of 10 MeV in the mw measurement. On the other hand, the fact that the individual 
sets predict uncertainties in the 10 MeV ballpark leaves hope that an improvement of the global 
PDF analysis will remove this bottleneck towards a precise mw measurement. 

The variation of the acceptance cut on the lepton pseudo-rapidity offers the possibility to 
scrutinize the dependence of the uncertainty on the flavor content of the proton and on the values 
of partonic-x. The preliminary results are not trivial, because of the correlations among the 
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densities enforced by the PDF sum rules. Increasing the value of the cut on \t]i\ reduces the PDF 
uncertainty on niw 
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